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1. Introduction 



There is considerable interest in the spectrum of D and D s mesons and 
of charmonium both theoretically and experimentally. 

On the theory side first principles calculations are usually lattice QCD 
fj . ■ computations (for recent work cf. e.g. [H El O HI O [6] ) . In the last couple of 

j — | years a lot of progress has been made, allowing the determination of hadron 

masses like the aforementioned mesons with rather high precision. For ex- 
ample 2+1 or even 2+1+1 flavors of dynamical quark are often used as well 
as small lattice spacings and improved discretizations, to keep discretiza- 
tion errors (in particular those, associated with the heavy charm quarks) 
under control. Some groups have even started to determine the resonance 
parameters of certain mesons from the spectrum of two-particle scattering 
$_i . states in finite spatial volumes (cf. e.g. |7J). 

Experimentally a large number of D, D s and charmonium states has 
been measured and additional and/or more precise results are expected in 
the near future both from existing facilities and facilities currently under 
construction, like the PANDA experiment at FAIR. Even though these ex- 
perimental results have been extremely helpful, to improve our understand- 
ing of QCD, they also brought up new and yet unanswered questions. For 
example the positive parity mesons Z)* and D s \ are unexpectedly light, 
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which is at the moment not satisfactorily understood and also quite often 
not reproduced by lattice QCD computations or model calculations. 

Moreover, performing a precise computation of certain meson masses is 
often the first step for many lattice projects not primarily concerned with 
spectroscopy. As an example one could mention the semileptonic decay 
of B and B* mesons into positive parity D mesons [8], whose masses and 
operator contents are an essential ingredient for any corresponding lattice 
computation. 

This is a status report about an ongoing lattice QCD project concerned 
with the computation of the spectrum of mesons with at least one charm 
valence quark. We present preliminary results for D mesons, for D s mesons 
and for charmonium states with total angular momentum J = 0, 1 and 
parity P = — , +. Parts of this work have already been published [9]. 

2. Simulation and analysis setup 

We use gauge link configurations with 2+1+1 dynamical quark flavors 
generated by the European Twisted Mass Collaboration (ETMC) [101 QTJ 
[T2l [TBI [T4"l I15| . Until now we have considered two ensembles (around 600 
gauge link configurations per ensemble) with (unphysically heavy) values for 
the light u/d quark mass corresponding to m n « 325 MeV, 457 MeV (lattice 
sizes (L/a) 3 x T/a = 32 3 x 64, 24 3 x 48). Our results are obtained at a single 
lattice spacing a ~ 0.086 fm. Consequently, a continuum extrapolation has 
not yet been performed. 

Meson masses are determined by computing and studying temporal cor- 
relation matrices of suitably chosen meson creation operators Oj. At the 
moment we exclusively consider quark antiquark operators. The quark and 
the antiquark are combined in spin space via 7 matrices and in color and po- 
sition space via gauge links (discretized covariant derivatives) such that the 
corresponding trial states Oj\Q) (|fi) denotes the vacuum) are gauge invari- 
ant and have defined total angular momentum and parity. Moreover, APE 
and Gaussian smearing is used, to optimize the overlap of the trial states 
Oj |0) to the low lying mesonic states of interest. More details regarding the 
construction of meson creation operators in twisted mass lattice QCD can 
be found e.g. in [16]. We plan to discuss these operators, their structure and 
their quantum numbers in detail in an upcoming publication. For the com- 
putation of the corresponding correlation matrices (O At) O(0)) we resort to 
the one-end trick (cf. e.g. p2]). Meson masses are then determined from 
plateaux values of corresponding effective masses, which we obtain by solv- 
ing generalized eigenvalue problems (cf. e.g. |18j). Disconnected diagrams 
appearing in charmonium correlators are currently neglected. 

For both the valence strange and charm quarks we use degenerate twisted 
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mass doublets, i.e. a different discretization as for the corresponding sea 
quarks. We do this, to avoid mixing of strange and charm quarks, which 
inevitably takes place in a unitary setup, and which is particularly problem- 
atic for hadrons containing charm quarks |14l [T5] . The degenerate valence 
doublets allow two realizations for strange as well as for charm quarks, ei- 
ther with a twisted mass term +ifi s , c j5 or —i(i s ,c75- For a quark antiquark 
meson creation operator the sign combinations (+, — ) and (— ,+) for the 
quark q and the antiquark q are related by symmetry, i.e. the correspond- 
ing correlators are identical. These correlators differ, however, from their 
counterparts with sign combinations (+, +) and (— , — ), due to different dis- 
cretization errors. In section [3] we will show for each computed meson mass 
both the (+, — ) = (— , +) and the (+, +) = (— , — ) result. The differences 
are 0(a 2 ) due to automatic 0(a) improvement inherent to the twisted mass 
formulation. These mass differences give a first impression regarding the 
magnitude of discretization errors at our currently used lattice spacing. 

Using (+, — ) = (— , +) correlators we have tuned the bare valence strange 
and charm quark masses /x s and \x c to reproduce the physical values of 
2m 2 K — m\ and mc, quantities, which strongly depend on fj, s and fj, c , but 
which are essentially independent of the light u/d quark mass. 

3. Numerical results 

In Fig. [JJwe present our results for the D and D s meson spectrum. For 
every state we show five data points: 
Red circles and crosses: 

lattice results at m n ~ 325 MeV, twisted mass sign combinations 
(+,-) = (-,+) and (+,+) = (-,-), respectively. 
Blue stars and boxes: 

lattice results at m^ « 457 MeV, twisted mass sign combinations 
(+,-) = (-,+) and (+,+) = (-,-), respectively. 
Gray triangles: 

experimental result from the PDG |19j . 
The differences between sign combinations (+, — ) = (— , +) and 
(+, +) = (—,—), which are <3%, indicate the magnitude of discretization 
errors at our currently used lattice spacing a « 0.086 fm. 

While for the negative parity states lattice and experimental results 
agree rather well, there is a clear discrepancy in particular for the positive 
parity D s states D* and D s \. Similar findings have been reported in other 
lattice studies, e.g. [TJE], and in phenomenological model calculations, e.g. 
|20j . This discrepancy might be an indication that these states are not 
predominantly qq states, but e.g. rather four quark states of molecular or 
tetraquark type. We plan to investigate this possibility within our setup in 
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Fig. 1. The D and D s meson spectrum. 



the near future. The necessary techniques have already been developed and 
recently been applied to light scalar mesons [21j . 

Another challenging, but important task is the separation of the two 
J = 1+ states, D 1 (2420), D x (2430) and D al (2460), D sl (2535), respectively. 
In the limit of a static charm quark one of these states has light cloud angular 
momentum j = 1/2, while the other has j = 3/2. To assign corresponding 
approximative j quantum numbers, when using charm quarks of finite mass, 
is e.g. important, when studying the decay of a B or B* meson to one of 
the positive parity D** mesons (which include the mentioned D\ (2420) and 
D\ (2430) states) in a fully dynamical setup (cf. e.g. [22], [23] for a recent lat- 
tice computation in the static limit). The correct identification of j = 1/2 
and j = 3/2 states can be achieved by studying the eigenvectors obtained 
during the analysis of correlation matrices; the largest eigenvector compo- 
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nents point out the dominating operators, which, after a Clebsch-Gordan 
decomposition into light and heavy angular momentum contributions, can 
be classified according to j = 1/2 or j = 3/2. 

In Fig. [2] we present our results for the charmonium spectrum. Be- 
cause of the two rather heavy valence quarks, we expect considerably larger 
discretization errors as for the corresponding D or D s meson states. The 
differences between lattice and experimental results are most prominent for 
the negative parity charmonium states (around 5%). We plan to explore in 
one of our next steps, whether discretization errors account for these dif- 
ferences by performing similar computations on ensembles with finer lattice 
spacings and by studying the continuum limit. 



> 

CD 



3.6 
3.5 

3.4 
3.3 
3.2 
3.1 
3 
2.9 
2.8 



charmonium 

XcO f 



Xc1 



jL 



m I =325MeV (+,-) ^-Q- 
m I =325MeV (+,+) ,— X- 
m I =457MeV (+,-) ^^ 
m^yMeV (+,+) ^4^- 
PDG > --A- 



Fig. 2. The charmonium spectrum. 
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